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The carbothermal processing of lunar regolith has been proposed as a means to produce 
carbon monoxide and ultimately oxygen to support human exploration of the moon. In 
this process, gaseous methane is pyrolyzed as it flows over the hot surface of a molten zone 
of lunar regolith and is converted to carbon and hydrogen. Carbon gets deposited on the 
surface of the melt, and mixes and reacts with the metal oxides in it to produce carbon 
monoxide that bubbles out of the melt. Carbon monoxide is further processed in other 
reactors downstream to ultimately produce oxygen. 

The amount of oxygen produced crucially depends on the amount of regolith that is molten. 
In this paper we develop a model of the heat transfer in carbothermal processing. Regolith 
in a suitable container is heated by a heat flux at its surface such as by continuously shining 
a beam of solar energy or a laser on it. The regolith on the surface absorbs the energy and 
its temperature rises until it attains the melting point. The energy from the heat flux is 
then used for the latent heat necessary to change phase from solid to liquid, after which 
the temperature continues to rise. Thus a small melt pool appears under the heated zone 
shortly after the heat flux is turned on. As time progresses, the pool absorbs more heat 
and supplies the energy required to melt more of the regolith, and the size of the molten 
zone increases. Ultimately, a steady-state is achieved when the heat flux absorbed by the 
melt is balanced by radiative losses from the surface. 

In this paper, we model the melting and the growth of the melt zone with time in a 
bed of regolith when a portion of its surface is subjected to a constant heat flux. The 
heat flux is assumed to impinge on a circular area. Our model is based on an axisym- 
metric three-dimensional variation of the temperature field in the domain. Heat transfer 
occurs only by conduction, and effects of convective heat transport are assumed negligible. 
Radiative heat loss from the surface of the melt and the regolith to the surroundings is 
permitted. We perform numerical computations to determine the shape and the mass of 
the melt at steady state and its time evolution. We first neglect the volume change upon 
melting, and subsequently perform calculations including it. Predictions from our model 
are compared to test data to determine the effective thermal conductivities of the regolith 
and the melt that are compatible with the data. 


Nomenclature 


C p Specific heat, J/(kg K) 

d n Normal displacement of the melt free surface, m 

h Convective heat transfer coefficient, W/(m 2 K) 

k Thermal conductivity, W/(m K) 

L Latent heat of melting, J/kg 

go Heat flux on the melt free surface, W/m 2 
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Q Rate of energy impinging on the melt surface, W 

r Radial coordinate, m 

Rbeam Radius of the beam of heat flux impinging the melt free surface, m 
iimeit Radius of the melt at its free surface, m 

t Time, s 

T Temperature, K 

T a Ambient temperature, K 

T m Melting point, K 

T s Temperature of the melt free surface, K 

v n Normal velocity of the melt free surface, m/s 

V Transformed temperature, K 

Knelt Melt volume, m 3 

Koid Volume of void space above melt, m 3 

X Location of the melting front, m 

Greek 

ct a bs Absorption coefficient of the melt surface 
p Density, kg/m 3 

e, e e ff Emissivity and effective emissivity of the melt and regolith surfaces 
n Thermal diffusivity, m 2 /s 

a Stefan-Boltzmann constant, a = 5.67 x 10 -8 W/(m 2 K 4 ) 

Subscript 

melt Pertaining to the melt 

reg Pertaining to the regolith 


I. Introduction 

A key element in NASA’s plans for exploration involves the production of oxygen for life support as well 
as use as a propellant from in-situ resources. For example, Sanders et al. 1 described NASA’s plans for in-situ 
resource utilization on the Moon. Various processes are being considered for the extraction of oxygen that is 
available in the lunar regolith as oxides of various materials. These processes include reduction using hydro- 
gen, molten salt electrolysis, and carbothermal processing. The carbothermal processing of lunar regolith 
using methane as the source for carbon was first performed by Rosenberg and co-workers at Aerojet in the 
1960s (Rosenberg et al.) 2 and developed further for lunar exploration by Orbital Technologies Corporation 
(ORBITEC). 3 

In the carbothermal process, gaseous methane is pyrolyzed above the hot surface of a molten zone of lunar 
regolith, depositing carbon on the surface and generating hydrogen which is carried away by the exiting 
stream of gas. The deposited carbon mixes with the melt and reacts with the metal oxides in it to produce 
carbon monoxide that bubbles out of the melt. Oxygen is ultimately produced by processing the carbon 
monoxide in reactors downstream. The amount of oxygen produced and its rate of production crucially 
depend on the volume of the melt. We have developed a chemical processing model of the carbothermal 
process 4 to predict the rate of production of carbon monoxide, and the volume of the melt is an input 
parameter to that model. 

In this paper we analyze the heat transfer aspects of the carbothermal process. Our goal is to predict 
the melt volume as it changes with time, as a function of the rate of heat input, surface radiative heat loss, 
and values of the thermophysical parameters, chiefly the thermal conductivity of the melt and regolith. We 
also determine the melt volume at steady state, if one exists. 
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II. Model 


A. Background 

Figure 1 shows a schematic of the carbothermal process. A model of the processing of lunar regolith by this 
method is described in Ref [4]. In this analysis we focus on the heating and melting of the regolith. The 
regolith is contained in a suitably sized reaction chamber. The surface of the regolith is subjected to a heat 
source such as from a laser beam or a concentrated beam of solar energy. In tests performed terrestrially, 
the power and beam size are in the hundreds of watts and a few centimeters in diameter, respectively. For 
space applications, the beam power is expected to be around a kilowatt. The beam size will be larger, with 
the heat flux expected to be comparable to terrestrial tests. In our model, we assume that the heat flux 
impinges on the regolith surface over a circular area and the beam intensity is uniform (or axisymmetric at 
least). Before the regolith is heated, we assume that its temperature is initially uniform. For a boundary 
condition on the chamber walls, we either assume an adiabatic condition or set the temperature to the initial 
temperature of the regolith. 

In order to minimize the radiative heat loss from the surface of the melt and the regolith, and maximize 
the volume of the melt for a given rate of heat input, a radiation shield is used that surrounds the top of 
the chamber. The radiation shield also provides an enclosure for containment of the gases (methane, carbon 
monoxide, etc.) involved in the carbothermal process. We do not include the radiation shield directly in the 
model. Its effect in reducing the radiative heat loss is incorporated by specifying an effective value (that is 
lower) for the emissivity of the surface of the melt and the regolith. While radiative heat loss is predominant, 
the model permits convective heat loss from the surface to the gas residing above it. 

Heat flux 



Figure 1. Sketch of the carbothermal process 

B. Governing equations 

We model the heating and melting of the regolith using the equations of heat transport by conduction and 
phase change. We neglect convective heat transfer in the melt. Convective flow in the melt is likely caused 
by temperature gradients due to natural convection and thermocapillary forces, as well as gravity-induced 
bubble motion when the reaction in the melt liberates gaseous carbon monoxide in it. Natural convection 
is likely to be weak since the melt is heated from above. When there is no reaction in the melt, bubble 
production is absent. Our model is restricted to such a scenario, and we will compare our predictions of 
melt mass with experiments without reactions. The effective themophysical properties of the regolith and 
the melt are treated to be different. A volume of regolith has a void fraction of nearly 40%. Thus the 
regolith properties are quite different from those of a bulk material of the same chemical composition. In 
our formulation of the governing equations, the regolith and the melt are treated as a single material with 
temperature dependent properties. The latent heat of melting is also accommodated by the single material 
formulation by an enthalpy method wherein the specific heat of the material is equivalently altered in a 
narrow temperature range around the melting point. 
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Since the void fraction in the regolith can be large, one can expect a large difference in the density be- 
tween the regolith and the melt. One complication due to the density change is the accompanying volume 
change upon phase change - the melt will occupy a smaller volume than the amount of regolith from which 
it is formed. Thus in addition to tracking the interface between the melt and the unmolten regolith (the 
melting front), the free surface between the melt and the gas phase above it must be tracked because of 
the volume change, which complicates the solution considerably. At first, we will analyze the heat transfer 
permitting the density change between the regolith and the melt, but without considering the volume change. 
We will include the effects of volume change later on. 

1. Time- dependent temperature field 

The energy equation is written as follows. 


pCp-^jr = V • (fcVT) (1) 

where p = p re g + (pmeit Preg ) t4(T T m ), k = k re g T (/c m eit k reg ) IA{T T m ), C p = C Preg T (Cp me it 

Cpre g ) U(T — T m ) + LS(T — T m ), and U(x) and 5(x) are the unit-step and delta functions, respectively, 
of argument x. In the computations, smooth numerical analogues for these functions are used. 

The initial condition is: 

T = T a at t = 0 (2) 

The boundary conditions are: 
dT 

= ^abs^o — e e ff cr(T 4 — T 4 ) — h(T — T a ) @ Regolith surface impinged by heat flux (3) 



— e e ff cr(T 4 — T 4 ) — h{T — T a ) @ Regolith surface outside heat flux 


dT 

-Tj— = 0 @ r = 0 (Axial symmetry) 


T = 


dT 

T a or — = 0 

on 


Container boundary 


( 4 ) 

( 5 ) 

(6) 


Once the solution for the temperature field is known, the melting front is determined as the surface where 
T = T m . The melt volume can then be obtained by a volume integral of the amount of material with 
temperature T > T m . 


2. Steady-state considerations 

We consider the question whether a steady state can be achieved where the melt volume and temperature 
field in the melt and regolith do not change with time. We assume that the container boundaries are 
insulated. For a steady state to exist, the total amount of heat input into the melt (which equals the heat 
flux impinging the melt that is absorbed minus the heat flux lost to the ambient from the melt surface by 
radiation and to a much smaller extent by convection) must be equal to the amount of heat loss from the 
non- molten regolith surface (due to radiation and convection). The energy equation at steady state can be 
written as 

V • (fcVT) = 0 (7) 

This equation is non-linear because the location of the melting front is unknown. Therefore, the regions in 
the domain where k = k me \ t and k = k reg are not known a priori. To overcome this issue, we use the well 
known Kirchhoff transformation. 


The energy equation is simplified to 


V(T) = 



V 2 y = 0 


(8) 

( 9 ) 
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The boundary conditions are 


= 0 @ Container boundary 


— &meit vr - = c^abs^o — e e ff cr(T 4 — T 4 ) — h(T — T a ) @ Regolith surface impinged by heat flux 
on 

d\^ 

—kme\t-Q— = — e e ff cr(T 4 — T 4 ) — h(T — T a ) @ Regolith surface outside heat flux 
where T in the above equations is the quantity 

r= bnelt y + 7; for _^L(T m -T a ) ( Qr T < Tm ) 

/v re g ^melt 

T = T m + V- p^(T m -T a ) for V > j^-(T m - T a ) (or T > T m ) 


While the problem is still non-linear because of the boundary conditions in Eqs (11) and (12), the energy 
equation has been rendered linear by the use of the Kirchhoff transformation. These equations are solved 
numerically for specified values of /q eg , &meit, T m , T a , <a a bs, <Tff, h and go- 


3. Volume change considerations 

As mentioned before, the regolith contains a substantial amount of void space and its density is quite a bit 
lower than a solid material of the same composition. The melt is expected to have a density similar to the 
solid material. Therefore, once the regolith melts, the volume occupied by the melt will be less than the 
volume of regolith from which it was formed. In the presence of gravity, the melt which is heavier sinks, 
and the melt free surface moves downward. The volume vacated by the heavier melt is occupied by the gas 
phase above it. 


Strictly speaking, in order to analyze the effect of the melt volume change, the equations of fluid motion 
need to to be analyzed along with the energy equation to describe the motion on the surface of the melt and 
within it. In our model, however, we use an approximate formulation without analyzing the melt motion 
in detail. We have used the conduction/phase change formulation of the energy equation (Equations 1 to 
6) where the location of the free surface between the melt and the gas phase above it is unknown and time 
dependent. In the numerical calculations, the free surface position is moved from its initial location by using 
either an interface velocity or displacement condition that are obtained by mass conservation upon phase 


change. 


(Pmelt Preg ) ^^melt 


v n = 0 for T < T r 


(/knelt Preg) 


dr, = 0 for T < T rr 


where V me \t and i? m eit are the instantaneous melt volume and melt radius on the free surface. These 
quantities are obtained by querying the instantaneous temperature field and demarcating the regions where 
the local temperature is greater than the melting point. As the free surface deforms the numerical grids 
are recalculated and the computations are carried forward in time using an Arbitrary Lagrangian Eulerian 
Scheme (ALE) with Winslow smoothing. 5 


III. Results 


The values of the thermophysical properties we have used in the numerical calculations are given in Table 
1. Many of the properties are not known precisely, and may not be the same for lunar regolith and its 
simulants such as JSC-1A. The properties used for the melt are the same as those for solid basalt rock from 
which regolith is formed. We have varied the values of the properties, chiefly the thermal conductivity, over 
a range in order to explore its influence on the melt mass. The lower value of the thermal conductivity of 
regolith is applicable when it is exposed to vacuum, such as on the lunar surface, and the value at room 
temperature and pressure for JSC-1A is around 0.2 W/(m K). 6 The computations have been performed for 
three-dimensional axisymmetric situations using the software package COMSOL 3.5a. Results are described 
below for the temperature distribution and melt mass at steady state, their transient evolution without 
considering volume change, and finally how the melt mass is influenced by volume change. We also compare 
the numerical results with specialized results for transient one-dimensional conduction with phase change. 
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A. Steady state 


Figure 2 shows the temperature distribution and the location of the melting front marked as a red line at 
steady state, when the container boundaries are adiabatic. Figure 3 shows similar results when the container 
walls are maintained at a temperature T = T a = 300 K. In each figure, we show the computed results with 
k reg = 0.01, 0.2 and 0.5 W/(m K), with fc me i t = 2.1 W/(m K) in all cases. The conductivity of regolith covers 
the range from the lunar surface exposed to vacuum to that in the presence of a gas at ambient temperature 
and pressure. Equations (9) to (14) using the Kirchhoff transformation have been used in these calculations. 

We see from Figures 2 and 3 that the thermal conductivity of regolith has a major role in the obtained 
melt volume, shown by the region with T > 1400 K (the melt masses obtained are given in the figure 
legends). For a fixed thermal conductivity of the melt, the melt volume is inversely proportional to the 
thermal conductivity of the regolith. At steady state, a fixed amount of heat that enters the melt must be 
lost from the cold regolith surface outside the area of the beam of incident heat flux. When the regolith 
thermal conductivity is low, the temperature gradient in it must increase to transfer the heat from the melt 
front (which is at the melting point) to the regolith surface (which is nearly at the ambient temperature). 
Thus the melting front location must be closer to the cold regolith surface when its conductivity is low, and 
therefore the melt volume is larger. It is also seen from Figures 2 and 3 that when the regolith conductivity 
is relatively large ( k reg = 0.2, 0.5 W/(m K)), the melt volume is less sensitive to the boundary condition on 
the container walls than for k reg = 0.01 W/(m K). 


0.05 


0.04 


0.03 


0.02 


0.01 


0 


-0.05 -0.04 -0.03 -0.02 -0.01 0 





(a) 


(b) 



( c ) 


Figure 2. Temperature distribution and location of the melting front at steady state. /c reg = 
(a) 0.01, (6) 0.2, (c) 0.5 W/{m K ), k me \ t = 2.1 W/(m K), a a b s = e e ff — 0.9, h = 0, heat flux beam diameter = 
0.75 in, Q=470 W, T a =300 K, container radius and depth = 5 cm, adiabatic container boundaries. The melt 
masses obtained are (a) 1066 (b) 55 and (c) 17 grams. 
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Density 

Specific heat 

Thermal 

conductivity 

Melting point 

Latent heat 


kq/m^3 

J/(kq K) 

W/ (m K) 

K 

J/kq 

Regolith 

1800 

840 

0.01 - 0.5 

1400 

435000 

Melt 

2850 

800 

1.5 - 2.1 


Table 1: Thermophysical properties of regolith and melt. 




(a) (b) 


□ 


-0.0S -0.04 -0.03 -0.02 -0.01 0 

Min: 300 

( c ) 

Figure 3. Temperature distribution and location of the melting front at steady state. /c reg = 
(a) 0.01, (b) 0.2, (c) 0.5 W/(m K ), k me \ t = 2.1 W/{m K ), o a bs — e eff — 0.9, h — 0, heat flux beam diameter = 
0.75 in, Q=470 W, T a =300 K, container radius and depth = 5 cm, container boundaries at T = T a . The melt 
masses obtained are (a) 474 (b) 31 and (c) 13 grams. 

B. Transient melt growth 

Figures 4 and 5 show the temperature distribution and the melting front location as they evolve with time 
for k reg = 0.05 and 0.5 W/(m K), respectively, with k me \ t = 2.1 W/(m K). We have used a somewhat larger 
value for k reg at the low end than in the results discussed at steady-state because for smaller values of k reg , 
the numerical calculations exhibited oscillatory temperature fields near the melting front, and the front shape 
was not smooth. At the final time step t = 4800 s, we see in Figure 4 that the melt growth is still evolving, 
while in Figure 5 for a larger regolith conductivity, the melt shape is quite close to a steady state. In both 
Figures 4 and 5, we see that melt radius increases with depth below the surface at first, attains a maximum 
radius, and then gradually decreases to zero at the location where the melt depth is a maximum. For both 
melts, these features of the melt shape are evident at all the times shown in the figures. The melt radius on 
the surface of the melt is constrained by the surface heat loss and is only slightly larger than the radius of 
the beam of heat flux that impinges the surface (beam diameter is 0.75 in). The surface melt radius achieves 
its steady value rather quickly (on the order of a minute) after the heat flux is turned on. 
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(C) (d) 


Figure 4. Temperature distribution and location of the melting front for a growing melt at (a) t=300 s, (b) 
t=600 s, (c) t=1200 s (d) t=4800 s. k re g = 0.05, k me \ t = 2.1 W/(m K), o a b s = e e ff — 0.9, h — 0, heat flux beam 
diameter = 0.75 in, Q=470 W, T a =300 K, container radius and depth = 5 cm, adiabatic container boundaries. 
A 3.5 x 3.5 cm region is shown. The melt masses obtained are (a) 20 (b) 31 (c) 46 and (d) 91 grams. 

Figure 6 shows the time evolution of the temperature on the surface of the melt at r = 0 for k reg = 0.05 and 
0.5 W/(m K). After the heat flux is turned on, the surface temperature rapidly rises from its initial value of 
300 K, and achieves the melting point of 1400 K in approximately 1 s. The surface temperature rises close 
to 2200 K in the time scale on the order of a minute, and thereafter changes slowly with time (T s = 2215 K 
at 100 s and 2238 K at 4800 s for k reg =0.05 W/(m K)). The maximum surface temperature that can ever 
be achieved is obtained by balancing the incident heat flux and the surface heat loss, ignoring conduction 
into the melt. Neglecting convective heat loss, the maximum surface temperature can be computed as 
T s = (<a a bs^o/( ecr ) T Tq ) = 2322 K. The ultimate surface temperature achieved (Figure 6) is somewhat 

less due to the heat flux conducted into the melt, not accounted for in the heat balance in computing T s . 

Figure 7 shows the predicted melt mass and how it compares with the masses obtained from test data. 7 These 
tests were performed in air and there were no chemical reactions in the melt. In obtaining the predictions, 
we have varied the thermal conductivity values of the regolith and melt in a small range to determine the 
sensitivity of the melt mass to these values. At t = 4800 s, the predicted melt mass is in the range from 
approximately 10 g (for k reg =0.5 and fc me i t =1.5 W/(m K)) to 30 g for k reg =0.2 and fc me i t =2.0 W/(m 
K)). It appears that for this test data, k reg =0.5 and fc me it =2.0 W/(m K) predicts the melt mass well for 
melt times greater than 1000 s. The value of k reg is somewhat larger than the measurements in Ref [6] at 
room temperature, but these authors have reported an increase in the apparent thermal conductivity due to 
stresses created by thermal expansion of the simulant particles against the rigid chamber wall. For reasons 
that are unknown, a lower melt conductivity (fc me it =1.5 W/(m K)) appears to fit the data better for earlier 
times. Perhaps k reg and fc me i t are temperature dependent, or the melts at earlier times are not well- formed 
and have regions of regolith particles trapped within them. In generating the results shown in Figure 7, we 
have assumed that the container walls are adiabatic and there is no convective heat loss from the surfaces 
of the melt and the regolith to the gas phase above them. We have used an emissivity of e = 0.9 for these 
surfaces, as a radiation shield was not used in the tests from which the data was obtained. 

C. Effect of volume change 

Figure 8 shows the time evolution of the temperature distribution, and the locations of the melting front and 
the free surface between the melt and the gas phase above it. The free surface motion is modeled using the 
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-0.025 -0.02 -0.015 -0.01 -0.005 0 -0.025 -0.02 -0.015 -0.01 -0.005 0 

Min: 323 Min: 499 


( c ) ( d ) 

Figure 5. Temperature distribution and location of the melting front for a growing melt at (a) t=300 s, (b) 
t=600 s, (c) t=1200 s (d) t=4800 s. k re g = 0.5, k me \ t = 2.1 W/(m K), o a b s = e e ff = 0.9, h — 0, heat flux beam 
diameter = 0.75 in, Q=470 W, T a =300 K, container radius and depth = 5 cm, adiabatic container boundaries. 
A 2.5 x 2.5 cm region is shown. The melt masses obtained are (a) 9 (b) 11 (c) 12 and (d) 14 grams. 


velocity boundary condition (Eq (15)). Similar results using a displacement condition for the free surface 
(Eq (16)) are shown in Figure 9. The predicted shapes of the void region at the rim of the melt region are 
different between these two results, with the velocity condition predicting a gradual shape change in the rim, 
and the displacement condition shows a sharper change. In both calculations we have used T — T a = 300 
K as the boundary condition on the container walls. Also, a convective heat transfer coefficient h m 10 
W / (m 2 K) has been used on the regolith and melt surfaces in contact with the gas phase. While the melt 
free surface moves, the heat flux is assumed to continuously shine on it and impinge perpendicularly with 
a constant beam diameter of 0.75 in (beam radius is 0.0095 m). Figures 8 and 9 show a large void volume 
(occupied by the gas phase) as the regolith melts, due to the change in volume on melting. Conservation of 
mass yields that (Kr 0 id + K ielt) Preg = Knelt /knelt => F void +F melt = p me it § = ~^850^ = 0-37. 

The shape of the predicted melt free surface in Figures 8 and 9 is worthy of discussion. In the presence of 
a gravitational body force, the free surface of a liquid is expected to be horizontal under static conditions. 
Deformations to the shape of the free surface can occur due to hydrodynamic stresses, liquid rotation, and 
capillary forces due to the effects of surface tension, especially near the meniscus region where the free 
surface comes into contact with a solid. Note that none of these effects have been included in our model. 
The melt free surface is not horizontal in the results shown in Figures 8 and 9 (the melt is the region with 
T > 1400 K) near its outer rim. While we have restricted the free surface deformation to be strictly vertical 
in using the conditions in Eq (15) and Eq (16), the results show that effect of gravity in maintaining a 
horizontal free surface is not strongly implemented by these conditions. The calculations shown in Figure 
9 with a deformation condition on the free surface seem to be better in this regard, and predict a smaller 
deviation from a flat shape for the melt free surface. We believe, in hindsight, that the approximate way in 
which we have used Eq (15) or Eq (16) to model the motion of the free surface to accommodate the volume 
change, without calculating the velocity field in the melt in detail (using the Navier-Stokes equations) and 
its effect on the heat transfer in the melt are responsible for the incorrect shape of the melt free surface that 
is predicted. 
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Figure 10 shows the predicted melt masses including the effect of volume change and how they compare with 
the test data. The velocity condition (Eq (15)) has been used in these results, but representative calculations 
show that the masses obtained using Eq (15) or Eq (16) are almost identical. Unlike in Figure 7, a single 
curve with fc reg = 0.5 and fc me i t = 2 W/ (m K) does not appear to represent the test data for t > 1000 s. 
This data is adequately bounded by predictions with fc me i t in the range 1.8 to 2.0 W/(m K). The data at 
earlier times is predicted accurately with fc reg = 0.5 and fc me i t = 1.2 W/ (m K), compared with fc reg = 0.5 
and fcmeit = 15 W/ (m K) in Figure 7. Note that in Figure 10 the container walls are at a fixed temperature 
and the surface convective heat transfer coefficient h = 10 W/(m 2 K), while in Figure 7 the container walls 
are adiabatic and h = 0. In general, the volume change appears to result in a slightly larger melt mass for 
a fixed fc reg and fcmeit compared to predictions that do not consider the volume change. Also, the predicted 
melt mass curves are flatter for larger times, indicating that the volume change tends to somewhat shorten 
the time taken to achieve steady state. It is likely due to the increase in radiative loss resulting from the 
increase in surface area of the melt and hot regolith surfaces that are evident in Figures 8 and 9. 

D. Comparison with one-dimensional calculations 

Carslaw and Jaeger 8 provide a solution for a one-dimensional phase-change heat transfer including the effects 
of conduction and convection induced by a density change. They assume that the temperature of the surface 
of the solid is held fixed at a constant value above the melting point, for all time. They derive an analytical 
solution for the location of the melting front as a function of time, which is given below. 
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The surface temperature T s is determined balancing the impinging heat flux and the surface radiative loss 
as 

1/4 

+X 4 ’ 
e a 


(19) 


rji | ^abs Q , /t-i4 

S “ U^beam " 

We find that the melt mass calculated from these equations is much larger (by about a factor of five at 
t = 4800 s) compared with the experimental data shown in Figure 7. This is likely due to the fact that the 
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Figure 7. Comparison of experimental and calculated melt mass 


one-dimensional model fails to capture the heat loss from the molten region to the regolith surrounding it. 
Furthermore, in the absence of the lateral heat loss, there is no steady state that is predicted. Therefore 
the one-dimensional analysis is of limited value and likely useful for only very short times after the start of 
melting. 

IV. Conclusion 

In this study, we have analyzed the time-dependent volume of a melt zone that is generated when a portion 
of the surface of a bed of regolith is subjected to a constant heat flux, and there is radiative heat loss from 
the melt and regolith surfaces. Under these conditions, we find a steady- state to exist with the melt volume 
being inversely related to the thermal conductivity of the regolith. The transient evolution of melting has 
been performed with and without considering the effect of volume change due to melting. The predicted 
melt masses have been compared to test data, and we find that thermal conductivities of 0.5 and 2 W/(m 
K) for the regolith and melt, respectively, adequately represent the data. 
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Figure 8. Temperature distribution and location of the melting front for a growing melt at (a) t=300 s, (b) 
t=600 s, (c) t=1200 s (d) t=4800 s. k reg = 0.5, fc me i t — 2.1 W/(m K), a a ^ s = e e ff = 0.9, h = 10 W/(m K), heat flux 
beam diameter = 0.75 in, Q=470 W, T a = 300 K, container radius and depth = 5 cm, container boundaries at 
T = T a . Velocity condition (Eq 15) used for the motion of the free surface. A 2.5 X 2.5 cm region is shown. 
The melt masses obtained are (a) 10 (b) 12 (c) 13 and (d) 14 grams. 


NAS A/TM— 20 1 2-2 1 7440 


12 






Figure 9. Temperature distribution and location of the melting front for a growing melt at (a) t=300 s, (b) 
t=600 s, (c) t=1200 s (d) t=4800 s. k reg = 0.5, /c me it = 2.1 W/(m K ), o a b s = e e ff = 0.9, h = 10 W / (m K), heat flux 
beam diameter = 0.75 in, Q=470 W, T a = 300 K, container radius and depth = 5 cm, container boundaries 
at T — T a . Displacement condition (Eq 16) used for the motion of the free surface. A 2.5 X 2.5 cm region is 
shown. The melt masses obtained are (a) 10 (b) 12 (c) 13 and (d) 14 grams. 



Figure 10. Comparison of experimental and calculated melt mass with the effect of volume change 
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